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Abstract 

This  paper  reveals  new  insights  in  the  aeroelasticity  and 
flight  mechanics  of  flexible  aircraft  by  obtaining  and 
solving  the  equations  of  motion  for  a flexible, 
accelerating,  rotating  aircraft.  We  illustrate  the  approach 
for  three  cases  of  increasing  complexity:  The  first  case 
is  a "sprung"  pendulum.  It  shows  when  rigid  body 
angular  velocities  can  be  important  in  the  flexibility 
equations  as  they  approach  as  the  flexible  frequencies. 
The  second  case  is  a typical  section  airfoil  on  an 
accelerating,  rotating  fuselage.  It  applies  Lagrange's 
equations  to  a longitudinal  problem  in  inertial 
coordinates,  then  transforms  the  equations  to 
noninertial,  body  - fixed  coordinates  for  solution.  It  also 
shows  when  rigid  body  rotations  and  longitudinal 
accelerations  must  be  included  in  the  flexibility 
equations.  The  third  case  is  the  general 
longitudinal/lateral  motion  of  an  accelerating,  rotating, 
flexible  vehicle.  Rather  than  setting  up  the  general 
problem  in  inertial  coordinates  and  then  transforming  to 
body  - fixed  coordinates,  instead  we  use  the  idea  of 
"quasi  - coordinates".  We  establish  a general  form  for 
Lagrange's  equations  in  the  noninertial,  body  - fixed 
coordinates.  The  paper  gives  the  general  equations  and 
reduces  them  to  a special  case  of  a "flat"  airplane.  It  also 
gives  guidelines  as  to  when  the  rigid  body  rotations  and 
accelerations  are  important  factors  in  the  flexibility 
equations. 

1.  Introduction 

For  many  years  there  has  been  a search  for  a practical 
set  of  "unified"  equations  of  motion  that  can  be  used  in 
all  of  the  disciplines  of  aerodynamics,  structures  and 
stability  and  control  of  flexible  aircraft.  Such  an 
approach  would  allow  the  customary  determination  of 
the  effects  of  structural  flexibility  on  aircraft 
performance,  stability  and  air  loads.  An  added  benefit  is 
that  it  would  also  allow  us  to  determine  the  effects  of 
the  "rigid  body"  motions  on  acroelastic  characteristics 
such  as  control  - effectiveness,  divergence  and  flutter. 
Further  (and  most  importantly),  it  would  allow  all  of 
those  engineering  problems  to  be  treated  by  subsets  of  a 
single  set  of  "unified"  equations.  In  effect  we  want  to 
convert  the  acroelastic  problems  into  coordinate  systems 
and  equations  that  arc  conventional  for  aircraft  flight 
mechanics,  stability  and  control. 


In  many  aircraft  applications,  the  mutual  coupling  of 
rigid  body  and  flexible  motions  has  been  small  because 
the  vehicle  angular  velocities  and  flexible  frequencies 
were  well  separated.  However,  there  have  been  recent 
examples  of  large  aircraft  where  flexible  frequencies  (say 
2 Hz)  begin  to  approach  the  rigid  body  angular 
velocities  (say  1 Hz).  Other  cases  have  been  known 
where  the  aerodynamic  forces  can  drive  the  structural 
frequencies  and  the  rigid  body  frequencies  close  together. 
In  both  cases  the  coupling  effects  should  be  accounted 
for  in  the  lowest  order  equations  of  motion  to  obtain  the 
correct  modeling. 

To  develop  the  necessary  equations  we  must  account  for 
the  fact  that  the  aircraft's  body  - fixed  coordinate  system 
is  not  (in  general)  an  inertial  system.  Dus  to  et  al  1*1, 
Bekir  et  al  l^.l  and  Waszak  and  Schmidt  f 31  are  a few 
examples  of  earlier  attempts  which  have  had  to  leave 
out  crucial  terms  or  were  difficult  to  implement.  This 
paper  shows  that  a practical  set  of  equations  for  general 
problems  is  available  through  the  use  of  energy 
methods,  Lagrange's  equations  and  "quasi  - coordinates". 

We  illustrate  the  approach  for  three  cases  of  increasing 
complexity:  The  first  case  is  a "sprung"  pendulum.  It 
shows  when  rigid  body  angular  velocities  can  be 
important  in  the  flexibility  equations  as  they  approach 
as  the  flexible  frequencies.  The  second  case  is  a typical 
section  airfoil  on  an  accelerating,  rotating  fuselage.  It 
applies  Lagrange's  equations  to  a longitudinal  problem 
in  inertial  coordinates,  then  transforms  the  equations  to 
noninertial,  body  - fixed  coordinates  for  solution.  It  also 
shows  when  rigid  body  rotations  and  longitudinal 
accelerations  must  be  included  in  the  flexibility 
equations.  The  third  case  is  the  general 
longitudinal/lateral  motion  of  an  accelerating,  rotating, 
flexible  vehicle.  Rather  than  setting  up  the  general 
problem  in  inertial  coordinates  and  then  transforming  to 
body  - fixed  coordinates,  instead  we  use  the  idea  of 
"quasi  - coordinates".  We  establish  a general  form  for 
Lagrange's  equations  in  the  noninertial,  body  - fixer! 
coordinates.  The  paper  gives  the  general  equations  and 
reduces  them  to  a special  case  of  a "flat"  airplane.  It  also 
gives  guidelines  as  to  when  the  rigid  body  rotations  and 
accelerations  are  important  factors  in  the  flexibility 
equations. 

The  equations  become  somewhat  more  complicated,  and 
it  is  useful  to  examine  them  in  three  stages.  First,  some 
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insight  is  available  via  a simplification  which  considers 
the  rigid  body  motions  merely  as  constant  parameters. 
There  the  rigid  body  motions  alter  the  flexible 
frequencies  of  vibration,  thereby  altering  aeroelastic 
stability.  Second,  a more  exact  approach  is  to  recognize 
that  the  flexibility  equations  have  some  of  the 
characteristics  of  Mathieu's  classical  ordinary  differential 
equation.  The  similarity  to  Mathieu's  equation 
introduces  the  possibility  that  the  coupled  rigid  - 
body/flexible  motions  can  be  unstable  within  narrow 
ranges  of  frequencies  and  amplitudes,  even  without 
aerodynamic  forces.  Third,  the  ultimate  procedure  is 
always  available  - the  simultaneous  solution  (perhaps 
numerically)  of  the  fully  coupled,  non  - linear,  rigid 
body  and  flexibility  equations  of  motion  in  body  - fixed 
coordinates. 

2.  Lagrange's  Equations 

If  the  inertial  coordinates  of  a dynamic  system  can  be 
represented  in  terms  of  N independent  generalized 
coordinates: 

X = X{qt)  E = Y{qi)  Z-Ziqj) 


Olsenl  1 showed  a related  (and  sometimes)  simpler 
approach,  noting  that  if  we  could  write  the  partial 
derivatives: 

y _ dX  y _ dY  y __  dZ 
* “ dqt  1 ~ dqt  i 3qt 

Y d2-*  y <>2Y  7 d2Z 

y dqfdqj  ij  dq^qj  V dqgdqj 

We  don't  need  the  often  tedious  expressions  for  the 
kinetic  energy,  and  the  equations  of  motion  take  the 
form: 

i miflj  + ii  ”kjrtj<ik  = Qi  +ft(: If)  - S; 

j=  i j=  \k  =i 

where: 

mij  = f(xixj  + Wj  + ZiZjjdm 
mass 

mijk  ~ f{^iXjk  + XiXjk  + XiXjk 
mass 


Lagrange's  equations  1 can  describe  the  motion  of  the 
system: 

d (JL\  dL  n 

^{dqj  dqt  Vi 

where: 

L = Lagrangian,  T - U 
Qi  = Generalized  Force 
T = Kinetic  Energy 
U = Potential  Energy 


For  simple  geometries,  it  usually  is  a straightforward 
matter  to  write  down  the  inertial  coordinates,  inertial 
velocities,  kinetic  and  potential  energies,  the 
Lagrangian,  and  the  various  derivatives.  For  complicated 
geometries,  the  process  can  become  tedious,  but 
Whittakcrl  ^ I showed  that,  if  the  kinetic  energy  can  be 
expressed  in  terms  of  the  coefficients  Wy 
n n 

i=  1 ,/T 


Then  the  equations  of  motion  can  be  written 


+ 

h i 


n n 

22\ 


Vkl 


'fjqk  = Qi  +i(f) 


EL 

eg. 


where  the  Christoff  el  symbol  is: 

\j  k]  = \ (dmi  dmH  dmh  \ 

[ i J dqk  dqj  dqt  ) 


Even  though  Whittaker's  and  Olsen's  expressions  look 
simple  in  principle,  in  practice  their  implementation 
can  be  qtiite  lengthy  for  complicated  geometries  with 
many  degrees  of  freedom.  The  development  of  the 
required  expressions  can  be  greatly  assisted  by  symbolic 
algebra  software. 

3.  Example  of  Coupled  Rigid  - Flexible 
Motions.  The  "Sprung"  Pendulum 

In  the  first  case  we  want  to  determine  when  "rigid  body" 
motions  can  have  important  effects  on  the  flexible 
motions.  Consider  the  "sprung"  pendulum  which  is  free 
to  rotate  or  oscillate  about  the  origin  in  the  x,  y (or 
r,6 ) plane,  but  which  also  am  tains  a radial  spring  of 
linear  stiffness  k (Figure  1).  We  will  refer  to  the  angular 
motion  as  the  "rigid  body"  motion  and  the  radial  motion 
as  the  "flexible"  motion.  Proceeding  through  the  usual 
processf  ^ 1 of  the  inertial  coordinates,  inertial  velocities, 
virtual  displacements,  kinetic  energy,  potential  energies 
(due  to  stiffness  and  gravity),  the  Lagrangian  is: 

L = ^m[r2  +r2d2)-^k{r-rkf  -mg{rsia6-Yg) 

From  I ^grange's  equations  the  radi  al  differential 
equations  is: 

r + («)o-02)r  = i + ^ -gsind 


We  also  can  obtain  the  angular  equation,  but  we  can 
always  interpret  it  as  the  angular  force  required  to 
produce  the  stipulated  motions. 

3.1  Rotation  at  Constant  Angular 
Velocity 
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In  the  first  case  we  stipulate  that  the  pendulum  moves 
through  a complete  circular  motion  at  a constant 
angular  velocity  of  to.  Thenfl  = (Ot , and  the  radial 
differential  equation  is: 

r + (coo  - ™2)r  = 75T  + £rk  ~ £sinR* 


Regardless  of  the  radial  force,  the  radial  response  acts  as 
if  the  natural  frequency  in  the  radial  direction  was 


effective 


1 


2 
w0 


to 


3.2  Simple  Harmonic  Rotation 


In  the  second  case  we  stipulate  that  the  pendulum 
oscillates  through  an  amplitude  Bq  with  a constant 

frequency  to . Then  0 = Bq  sin  (Ot , and  the  radial 
differential  equation  is: 


, r 2 

ir+jto^ 


2a2 
a)  6i 


— (l  + cos2tof)]r  = 


+ 


The  complete  solution  of  the  radial  equation  depends  on 
the  LHS,  RHS  and  initial  conditions.  The  LHS  can  be 
converted,  with  a change  of  variables, 

, 1 * 1 /j2 

r = art  h = — --jBq  o = ifOo 

to  the  classical  Mathieu's  equation. 

r"  + (a-  2fccos2r)/-  = 0 


Mathieu's  equation  applies  to  the  vibrations  of  spinning 
satellites,  buckling  of  beams  with  periodic  end  forces, 
the  saturation  of  loudspeakers,  tides  in  circular  bodies  of 
water  and  many  other  problems.  In  our  application,  if 
the  radial  force  does  not  depend  on  r then  the  stability  of 

. two 

the  solutions  depends  only  on  the  frequency  ratio  — 
and  the  angular  amplitude  Bq 


Intuitively,  one  would  expect  that  the  effects  on 
stability  would  be  small  unless  the  angular  amplitude  is 
large  or  the  frequency  ratio  is  near  l.  Figure  2 (from 
McLachlan^)  shows  the  classical  plot  of  the  regions 
of  stability/instabilily  for  periodic  solutions  of 
Mathieu's  equation.  Regions  of  instability  are  shown  to 
be  emanating  from  the  points 


= 1>22,32, 


2 


n 


So  a question  becomes  - what  practical  values  of  a,  b 
put  the  solutions  into  the  stable  or  unstable  regions. 
For  instance,  in  the  neighborhood  of  a=l,  we  can  use 
Mcluchlan'sf^l  boundaries  to  obtain  a region  for 
instability  for  small  b.  Figure  3 shows  the  lower  and 


upper  bounds  of  the  narrow  unstable  region  for 
frequency  ratios  near  1 and  for  angular  amplitudes  up  to 
10  degrees.  The  instability  range  continues  to  widen  for 
higher  values  of  the  angular  amplitude.  For  frequency 
ratios  near  2, 3, 4,...,  the  instability  ranges  exist,  but 
over  ever  narrower  ranges  of  angular  amplitude. 

We  also  can  integrate  the  equation  numerically.  Figure 
4 shows  typical  time  histories  for  a frequency  ratio  of 
1.1,  damping  of  0.02  and  amplitudes  of  0.53, 0.54  and 
0.55  radians.  The  solutions  are  stable  for  0.53  radians 
and  unstable  for  0.54  and  0.55  radians.  Figure  5 gives  a 
general  pattern  for  the  smallest  amplitudes  to  produce 
instability  for  damping  of  0.02  and  frequency  ratios  up 
to  2. 

In  summary,  the  problem  of  the  "sprung"  pendulum 
shows  that  rigid  body  motions  can  affect  the  flexible 
motions: 

a.  Constant  angular  velocity  reduces  the 
"effective"  radial  natural  frequency; 

b.  Forced  sinusoidal  angular  motion  can 
produce  radial  instability  near  integer  values  of  the 
frequency  ratio  as  the  angular  amplitudes  grow  large. 

4.  Typical  section  Airfoil  on  an 
Accelerating.  Rotating  Fuselage 

The  second  case  is  a problem  that  is  closer  to  practical 
interest  - a typical  section  airfoil  on  an  accelerating, 
rotating  fuselage.  We  will  apply  Lagrange's  equations  in 
inertial  coordinates,  then  transform  the  equations  to 
noninertial,  body  - fixed  coordinates  for  solution.  We 
want  to  show  when  rigid  body  rotations  and 
longitudinal  accelerations  must  be  included  in  the 
flexibility  equations. 

Consider  a slender  airfoil  which  is  mounted  on  a slender 
fuselage.  (Figure  6).  The  fuselage  has  inertial 
coordinates  Xq  = Cj\,  1q  - (J2  ;uKl  pitch  angle 
B = q3.  The  airfoil  is  located  at  fuselage  position 
X = Xw  and  has  its  own  degrees  of  freedom  in  vertical 
translation  h = q4  and  rotation  6 = q$. 

4.1  Equations  of  Motion 

For  a general  point  on  the  slender  fuselage  and  airfoil 
the  inertial  coordinates  are: 

+ XC3  Y=q2+xs3 

+ xwc3-q4s3  + %c35 
+ xws3  + qAc3  +§s-35 


Fuselage  X = q4 

X = q\ 

Airfoil 

y-q2 
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Proceeding  through  the  usual  process  of  the  inertial 
coordinates,  inertial  velocities,  virtual  displacements, 
kinetic  and  potential  energies  (due  to  stiffness  and 
gravity)  and  the  Lagrangian  for  the  complete 
fusclage/airfoil  system,  we  establish  the  complete 
nonlinear  equations  in  inertial  coordinates.  Then,  it  is 
more  convenient  to  actually  solve  the  problem  in  a 
body  - fixed  coordinate  system.  So  we  define  the 
"apparent"  body  - fixed  components  of  the  vehicle's 
velocities  and  accelerations  by 

9*  = W3  + <72*3  qy  = 02<3  ~ 9l*3 


clx  "=  <?lc3  + Qls3  Qy  = <12C3  ~ 9bv3 


The  total  solution  requires  the  solution  of  the  five 
coupled  equations.  It  can  be  convenient  to  separate  the 
complete  equations  into  separate  rigid  equations  and 
flexible  equations.  If  we  separate  the  purely  rigid  body 
terms,  coupling  terms  and  purely  flexible  terms  - the 
linearized  flexibility  equations,  in  body  - fixed 
coordinates,  (with  damping  added)  are : 

\mw  \c4 

L^W  Av-Jl©] + [°  c5jl?5j 


t\k4  01  .2f/n* 

TLo  *s|  93 [sw 


Sw  1 J° 

tjJi* 


[04  1 

f mW  mWXW  ■*"  1 
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[ ^ W + XV,  SW  \ 

h 

— 8C '3 

Sw 

r Lo  AsjUJ 

We  simplify  further  by  dropping  second  order  terms  in 
<74,  <75  to  obtain  the  linearized  equations  in  terms  of 
noninertial.  body  - fixed  coordinates. 
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where 

mw  = Jdm 
wing 

Sw  = pdm 

wing 

lw  = f%2dm 

wing 

M33  = lf  + lw  + ™w*w  + 2SWXW 
A/34  = mwxw  “*■  Sw 

M35  = I W + Swxw 

M45  = sw 

M334  = mw04  + Swdfi 
^335  = Sjc/4  ~ xwQ 5) 


MA35  = Swd5 

My  3 = Sy  + xwmw  + Sw 


4.2  Separate  Rigid  and  Flexible 
Equations  - Body  Coordinates 


4.3  Vibration  Solutions 

Immediately  we  can  sec  hints  of  the  effects  of  the  pitch 

rate  ^3  and  the  acceleration  along  the  body  axis  qx  as 
they  alter  the  "effective  stiffness"  in  the  flexibility 

equation.  Assuming  <73,  qx  are  constants,  figures  7 and 
8 show  the  effects  of  aircraft  pitch  rate  on  the  coupled 
(unbalance  not  equal  zero)  airfoil  frequencies  of 
translation  and  rotation  for  the  airfoil  slightly  aft  and 
slight  forward  (1  chord)  of  the  aircraft  axis. 


In  the  case  of  the  unbalance  equal  zero,  the  vibration 
equations  are  uncoupled  and  the  translation  mode  just 
acts  as  if  the  "effective"  stiffness  is 

k\f  — k4~  mwq2  = k4  - (-g-)  j 

In  the  case  where  the  93  ~ ^ but  the  unbalance  is  not 
zero,  the  equations  remain  coupled,  but  the  torsion 
equation  acts  with  an  "effective"  stiffness 

h-S^-kJ 

V rw  cw  5 « / 


Whether  the  effective  torsional  stiffness  is  slightly 
larger  or  smaller  depends  on  the  sign  of  the  unbalance 
and  whether  the  aircraft  is  accelerating  or  decelerating. 


4.4  Aeroelastic  (Hypersonic! 
Equations 


We  use  (for  convenience)  hypersonic  aerodynamics  from 
piston  theory^  \ to  obtain  the  hypersonic  flexibility 
equations. 


sjHn  LiA«s«  ° 1 [94] 

^JK)  1 1 0 *3*5%]  ’Hd  AjJ l^s J 


JM  0 1 jjJO  I|  2 1 1 

J1  “(x-  + 6)  1 


*.  1 . I J°  °1U?4| 


.I1  lr«i_r  1 4+i,  ipn_.  r 1 1 (4  0 lfivj 

[h  -{XV^  + > [*w  3|  Av  J [0  j ~r  i 


where  the  nondimensional  variables  are  defined  by: 
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lengths 

c=chord  l=span 

* 

positions/coordinates  Xw  = CXW 

<?4  = ^4  = ^?jc  Qy  = cQy 

inertias 

c , 2-2 
stiffnesses 

2 2 2~2  2 

&4  = ^,04  /C5  = /wc»5  = mwc  rw(»s 

dampings 

c4  = 2^0)4^  c5  = 2mM,c2f^fi>5g'5 


time  and  frequency 


Cfio 


1 C(U 
A — y 


time  derivative 


Therefore  those  rigid  body  motions  also  can  alter  the 
aircraft  speeds  for  aeroelastic  divergence  and  flutter.  For 
example,  under  the  assumed  conditions  we  can  obtain  an 
approximate  expression  for  the  [A  required  for 
aeroelastic  divergence: 


which  shows  the  importance  on  the  divergence  speed  of: 

a.  £ , the  ratio  of  the  pitch  rate  to  the 
uncoupled  translation  frequency; 

q” 

b.  , the  relationship  of  the  acceleration  to 

the  torsional  frequency  and  the  radius  of  gyration 

c.  (j) , the  ratio  of  the  uncoupled  translation  and 
rotation  frequencies 

d.  , the  relationship  of  the  unbalance  to  the 
radius  of  gyration. 


air  density  and  gravity 


pcJc 


We  also  can  use  the  hypersonic  aeroelastic  equation  to 
do  an  eigenvalue  calculation  (dropping  the  RHS)  to 
obtain  flutter  solutions.  Figures  9 and  10  show 
representative  effects  of  pitch  rate  and  acceleration  on 
hypersonic  divergence  and  flutter  boundaries. 


aerodynamic  coefficients 

Qa=Nw  Qs  = Mw  = pV^C  CMw 

Geometry  integrals 

fl-ite-0.5  ij-fg-i  +0.333... 

As  usual,  the  damping  is  modified  with  the  aerodynamic 
damping  and  the  stiffness  is  modified  with  the 
aerodynamic  stiffness.  Flowever,  the  stiffness  also  has 
terms  that  are  proportional  to  the  nondimensional  pitch 

rate  <?3  and  the  nondimensional  acceleration  Qx  ■ 

4.5  Rigid  Body  Motions  as  Constant 
Parameters  in  the  Aeroelastic  Equations 

Now  again  consider  the  pitch  rate  and  the  aircraft 
acceleration  as  constant  parameters.  From  the 
differential  equation  and  the  vibration  solutions,  we 
know  that  the  rigid  body  pitch  rate  will  decrease  the 
bending  frequency  (even  if  only  slightly),  and  that  the 
rigid  body  acceleration  (or  deceleration)  along  the  body 
axis  ran  increase  or  decrease  the  torsional  frequency. 


4.6  Forced  Rigid  Body  Motions  in  the 
Aeroelastic  Equations 

Rather  than  assume  that  the  rigid  body  motions  are 
constant  parameters,  we  can  assume  representative 
forms  for  their  time  dependent  motions  and  then  plug 
them  into  the  flexible  equations  of  motion.  We  need  the 

terms  #3  > % > Qx’Qy  an£^  *?3-  Following 

191 

Etkin’s1  1 notation  we  can  assume  the  time  dependent 
forms  for  the  oscillatory,  damped  speed,  pitch  angle  and 
angle  of  attack,  wherein  each  expression  the  terms  a and 
b are  assumed  constants: 

Q x = M0  + E(au^l  + buCij 
0 = #3  = 00  + E{adSl  + bpCi) 
a = Ot()  + + ^aQ) 

where 

Sn  = sin  nmrf,t  Cn  = cos  nm^t  E =e~  0>rh^rht 
(0rl,  = assumed  " rigid  body"  frequency 
grl,  = assumed  " rigid  body"  damping 

Noting  that  Qy  “ Qxa  and  combining  angles  where 
possible,  we  obtain  the  "forcing  terms"  to  be  included 
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in  the  equations  of  motion.  The  terms  cl3  ■>  0 x appear 
on  the  LHS  of  the  equation  and  influence  stability: 


The  major  point  is  that  terms  s->  <=  sin  XikiJ  ;ind 

Ci  = cos  2 gjrbt  on  the  LI  IS  introduce  behavior  like  that 

of  Mathieu's  equation. 

5.  The  General  Case  - Three  Dimensional 
Motion  of  a Flexible  Vehicle 

5.1  Geometry 

We  start  with  an  inertial  X,  Y,  Z coordinate  system  and 
a noninertial  X,  y,  Z system  that  can  accelerate  and 
rotate  in  the  X,  Y,Z  system  (Figure  1 1).  The  origin  of 
the  X,  y,  z is  located  in  the  inertial  system  at: 

X = X0  = qi  F = F0=%  Z = Zq  = <73 

The  orientation  of  the  X,  y,  z.  system  is  given  by  the 
conventional  sequence  of  Euler  rotations: 
ip  =g4  6 =q5  (j)  = q6 

5.2  Inertial  Coordinates 

Then  inertial  coordinates  of  a general  point  in  x,  y,  z 
arc: 

{X,Y,Z}  = {0i,?2.?3}  +[>]{*,w4 

where  [t]  the  Euler  transformation^  is  the  product 
of  three  transformations  that  depend  on  the  Euler  angles: 

[r]=[r4lr5]r6] 

\C4  -54  01  fc5  0 *1  ft  0 0 1 

[t4]-\m  c4  °l  ['%]_l  0 1 0l  [tg ] ==  1 0 c6  -i6l 

[0  0 ij  l-s5  0 c5J  [0  *6  c6J 

s;  = sin  q,-  c = cosqf 

We  also  write  the  "local"  coordinates  in  terms  of 
additional  generalized  coordinates  qj,q#...qn: 
n 

{x,y,z}  = J \{xi,yi,Zi}qi(t ) 
i=7 

to  obtain  the  inertial  coordinates  in  terms  of  the 
generalized  coordinates: 

n 

{X,  Y,  Z}  = {0!  ,02.9?}  + {x'  ’ Yi'zi  h ( ^ 

r=7 


The  kinetic  energy,  in  terms  of  the  inertkd  coordinates 
is  a lengthy  expression  which  shows  why  it  can  be 
useful  to  use  Olsen’sl  ® 1 form  of  the  equations  of 
motion  (which  doesn't  require  the  kinetic  energy),  rather 
than  Whittaker's^  I form  which  does  require  the  kinetic 
energy. 

5.3.  Overcoming  the  Tedious  Aspects 
- Quasi  Coordinates 

Again,  we  could  use  Eagrange's  equations  on  the 
I agrangian  in  inertial  coordinates  to  obtain  the 
equations  of  motion  for  the  flexible  system.  We  would 
be  accurately  accounting  for  all  of  the  inertia  couplings 
that  arise  from  the  fact  that  the  noninertial  x,  y,  z 
system  is  accelerating  and  rotating  in  the  inertial  X,  Y, 
Z system.  We  could  solve  the  problem  in  terms  of  the 
inertial  translations  q\ , t/2  ’ mk1  the  Euler  angles 
q4 , q^ , q^  and  then  transform  the  results  to  the 
translations  along  the  body  axes  qx,  qy,  qz  and  the 
instantaneous  angular  velocities  < ox,(Oy,mz , using 
the  transformations: 

{9jr>9y.9z}  =[ A1 

^wx,(Oy,  Mz  I = [ tr]  jc/4 , <75,  i/(, } 

where 

r -*5  0 1 1 

[a]=lc5A’6  c6  ol 

[C5C6  ~s6  0j 

'Ihe  approach  is  correct  in  principal.  However,  it  works 
easiest  for  special  cases  like  rotation  about  one  axis 
(where  the  time  derivative  of  the  appropriate  Euler  angle 
is  indeed  the  angular  velocity).  However,  it  suffers  from 
two  shortcomings  in  the  general  case  of  three 
dimensional  motions. 

First,  the  generalized  coordinates  t/j , > <73 

are  the  translations  in  the  directions  of  the  inertial 
coordinates.  We  would  like  to  replace  them  with  the 
translations  in  directions  of  the  noninertial,  body  - axis 
coordinates  qx,  qy,  qv 

Second,  the  generalized  coordinates  q4 , C/5 , 95 
are  the  Euler  angles.  Their  time  derivatives  q4 , <75 , t/5 
may  not  be  the  physical  angular  velocities  of  the  x,  y,  z 
system  for  general  motions.  We  would  like  to  replace 
them  with  the  physical  angular  velocities  of  the 
noninertial,  body  - axis  coordinates,  (Ox,  (Oy,  (O z. 
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However,  a much  more  elegant  and  simple  method  is 
available,  the  method  of  quasi  - coordinates  due  to 
Hamel  t * ® I and  Boltzmann  1 1 * ! Hie  term  "quasi  - 
coordinates"  refers  to  the  fact  that  we  cannot  (in  the 
general  case  of  three  dimensional  motions)  direedy 
integrate  the  angular  velocities  to  get  the  generalized 
coordinates.  Actually  all  we  are  doing  is  performing  the 
transformations  before  we  apply  Lagrange's  equations  to 
obtain  the  differential  equations,  rather  than  after  we  get 
them. 

Whittaker!  ® 1 and  Meirovitchf  1 ^ ] explain  the  method 
of  quasi  - coordinates  for  the  special  case  of  rotational 
motions.  Several  others,  among  them  Nayfeh  and 
Mookl  13],  give  applications. 


f ah  1 

uli 

\ at  1 

| | 

I <3?4  1 

| d«>x  I 

(04  } 

d j dL  1 
3?  ' {fat  ■ 

-\PUQs\ 

\j£\ 

[ *»z  J 

1 dL  1 

\jL_\ 
[<?®z  J 

[ft] 

'Hie  basic  idea  is  that  we  want  to  write  Lagrange's 
equations  in  a form  that  treats  directly  the  body  axis 
translations  qx,  qy,  qz  and  the  true  angular  velocities 

cox,  a)y,(Oz.  We  start  with  the  usual  form  of 

Lagrange's  equations  in  terms  of  the  original, 
independent  generalized  coordinates  q^ , <g>  - - ■ qn  and 

their  time  derivatives  4b  92  • ■ ■ ’ 


(JL\_  1 

(JL\ 

\*4i)  1 

{ dq  j 

Flexible  DOFs  remain  the  same  (except  that  we  must 
use  the  modified  I ngrangian  L ) 


d dL 

3F  dq. 


dL 

dQi 


Qi  for  i a 7 


where: 


r o 


-wz  a )y  1 

0 -wx  I 

(Dx  0 


The  Ingrangian  can  be  written  in  the  usual  form  in  the 
original  inertial  coordinates: 


L = 


9l>92>93> 94 ’95’ 96 ’97 v "Qn >\ 
V9l’92’93’  94’95’96’97’-  "Qn  / 


If  we  note  that: 


where 


{94.95.%} 

'y,  ^ 

(o  c$\s6 

C5C6 

M = [«]_1  =0  c6 

~s6 

1 *5*6 

*5C6 

we  can  obtain  the  equivalent  form  of  the  Lagrangian  in 
the  quasi  - coordinates: 

~( Qx’Qy,Qv’  94’  95’96’97v  "Qn ’ \ 

L = L| 

\Qx'Qy->Qz'(0X'(ay->(0z'Q7','"Qn, 

Then,  following  Whittaker!  ^ 1 we  can  obtain  the 
equations  for: 

Translation  DOFs: 


These  are  the  equations  of  motion  in  terms 
of  quasi  - coordinates.  They  are  the 
fundamental  advance  which  allows  us  to 
formulate  a unified  set  of  equations  that  can 
be  used  without  simplification  for  the 
aerodynamics,  structures  and  stability  and 
control  of  flexible  aircraft  - they  allow  us  to 
place  the  aeroelastic  problem  into  a 
coordinate  system  and  notation  that  is  used 
in  flight  mechanics  and  stability  and  control. 

5.4  Energies  - Noninertial  Body  - 
Axis  Coordinates 

The  kinetic  energy  in  terms  of  the  body  - axis  variables 

is: 

+ Sx[qya)z  - qzmy)  + Sy(qzmx  - qxa>z)+ Sz[qxmy  - qymx ) 
+‘rIxx[wy  ■*" ®z  ) — Ly  ®x IJ> y — f xz®x®z 
• +jlw(ya>l  +a>l ) -IyZwycoz  + (a£  + ay  ) 

+ + Syi/y  + Szqz 

+{*xy  “ 4y)®z  + i/yz  ~ lyz)mx  + Ozx  ~ ^ix)wy 
"*“2  (fix  + fyy  +fzz  ) 
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where  typical  inertia  integrals  are: 

IS 

II 

H 

velocity  along  the  x axis 

M-j  dm 

ii 

< 

velocity  along  the  y axis 

body 

4x  = w 

velocity  along  the  z axis 

Sx=Jxdm  = ^[J xidm)pi  = J Sx.  qt 
i i 

0)x  - P 

angular  velocity  about  the  x axis 

Sx=Jxdm  = '2  {jxidm)pi  = J Sx.  % 

My  ~ Q 

angular  velocity  about  the  y axis 

i i 

hy  = f xydm  = ^ (/  xiyJdm)cliclj  =2/x,>j  q^J 

body  ‘J  U 

(»z  = R 

angular  velocity  about  the  z axis 

ix,y,t)^^zi(x,y)qi(t) 

i 

hy  = f xydm  = ^(f '■*iyjdrn}ii4j  = E V,  qiq-> 
body 


id 


hj 


The  potential  energy  due  to  gravity  will  come  from  our 
gravitational  model.  In  the  case  of  a "flat  earth": 

Vg  = ~g  f(z-  Zref)dm 


, ^ ■ N 

--dy<Ii-Zref)m-  g[0  0 lj[r 


We  expect  the  potential  energy  due  to  flexibility  to  be 
of  the  form 

i=  Ij  =7 

Or  perhaps  a more  general  expression  for  larger 
deflections 

vbifl$j + ifii  vbijmfli 

i-1  j=l  i=l  j=7k~l 

The  author  has  performed  those  operations,  and  the 
complete  set  of  differential  equations  is  available  (but 
too  lengthy  to  present  here). 


The  equations  of  motion  become 


Translation 

m [°  0 -*]n  ,JQ  i 

Miv  |+l  o o sx  -wraf  + \ pq  -{p1  + 

\w\  [sy  -Sx  0 ][*]  \vP-UQ]  | PR  QR  jl  5 

t 55 1 Jn  m f (?.♦*)  1',  km 

M '(M  l o j [>2+o2)jJ  ‘ UJ 


Rotation 

i c o r ^ -/»  c ][/>i 

i c o -s,h  rM-v  c Nfil 

[-Sj  Si  0 J[Wj  | 0 C 

r 0 » G*lfi«l  \VQiWR  -OQ  1 j t C -cjttl 

il-P*  -QR  C l</»4<l  VP  unm\  ' +S!c«  0 I ' 

I PQ  fl3-'’2  -eeJl^J  I -w  -M2  JI>J  [-c3«  “ *5  J 

ffO  C J r-v  0 c , -(*+  PQ]  [#-<?','  1(6,  1 

" 0 0 I4/+V-W)  (p2-*2)  -(*-Pe)«,js,  \ 

i||0  C Oj  |C  -2 P -7Q\  j(W+ VG)  -(P-fl*)  -((2+  MT)j  J [S,  J 

\<SH  o d f:S',  ] [2«i)><4/(i*-!3Kj|  [Afj] 

4»l«  C Ol^S,  rti-W,') 

(C  0 0j'[Sj,|  “[  0 j [*)J 


Flexible 

(W+VP  - UQfc.  - (Q  - PR)Sxl . +(/'  + e*)5,( 

+2Mp2  -Qi 

j 


5.6  Simplification  to  the  "Flat" 
Airplane 

If  we  specialize  the  general  body  to  consider  an 
essentially  "flat"  surface  in  the  xy  plane  (Right  handed 
xyz  coordinate  system  attached  to  the  body),  Figure  12, 
and  make  the  usual  definitions 


The  equations  above  are  the  equations  to  solve  for  the 
static  and  dynamic  response  and  stability  of  a flexible, 
"flat"  aircraft  under  steady  state  flight  or  in  accelerations 
and  maneuvers.  They  are  nonlinear  and  mutually  couple 
the  overall  rigid  body  motions  with  the  flexible 
deflections.  They  can  be  used  for  analyses  of  aircraft 
performance,  stability  and  control,  flight  loads,  control 
effectiveness  and  aeroelastic  divergence  and  flutter.  Of 
course  they  are  more  complicated  than  the  conventional 
nonlinear  equations  for  rigid  body  motions  or  the  linear 
equations  for  aeroelastic  response  (which  are  coupled  to 
the  rigid  body  equations  only  through  the 
aerodynamics). 


5.7  Perfect  Masses  and  Modes 


In  the  special  case  of 

a.  Mass  symmetry  about  the  y axis 
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b.  Origin  at  the  center  of  mass 

c.  Perfectly  orthogonal  free  - free  modes 

the  equations  of  motion  simplify  further  to  the  usual 
rigid  body  equations  and  the  modified  flexible  equations: 


Translation 

[t/|  fWfi-VRl  r s$  1 fj 

\ V !-  + \UR- WP  i+  g\ -c5s6  h jf  \Fy  !• 
[w\  [VP-UQ\  l-c5c6J  [fJ 


Rotation 


1 p 

f QR 

\Q 

l i 

+ 

l*J 

yy  r</-\ 

I +1 

xx  yy 

where 


Flexible 


= Qi 


a\  = natural  frequency  of  the  perfect  ith  mode 


If  we  assume,  for  the  moment,  that  the  angular 
velocities  P and  Q are  constant,  then  one 
approximation  would  be  to  treat  the  flexible  equations 
as  if  the  effective  structural  frequency  for  any  mode  is 
just  replaced  by 


(i^+e2) 


On  the  other  hand,  since  P and  Q will  be  functions  of 
time,  the  actual  behavior  will  be  more  like  the  behavior 
of  solutions  to  Mathieu's  equations. 

Many  recent  developments  numerically  integrate  the 
linearized  equations  of  motion  with  nonlinear 
aerodynamics  on  the  RHS.  It  seems  that,  once  the 
analyst  has  committed  to  numerical  integration  of  the 
equations  of  motion,  there  is  very  little  additional  labor 
(or  computational  time)  to  use  the  more  comprehensive 
equations  of  motion  above. 

6.  Summary.  Conclusions 


coupling  mechanism  between  the  rigid  body  and  flexible 
motions,  even  in  the  absence  of  aerodynamics.  It 
appears  that  if  a flexible  frequency  is  up  to  1.3  - 1.5 
times  a rigid  body  frequency,  then  those  coupling  effects 
should  be  considered.  In  some  cases  ("slender"  aircraft) 
the  natural  frequencies  may  already  be  in  those  ranges. 

In  other  cases  (the  X - 29)  the  aerodynamic  forces  drive 
some  of  the  flexible  frequencies  down  toward  the  rigid 
body  frequencies. 

3.  The  airfoil  on  an  accelerating/rotating  fuselage  shows 
that  the  effective  bending  stiffness  is  reduced  by  a 
constant  pitch  rate.  It  also  shows  that  torsional  stiffness 
is  increased  or  decreased  by  constant 
acceleration/deceleration,  depending  on  the  sign  of  the 
unbalance.  The  results  modify  the  divergence  and  flutter 
speeds.  If  we  impose  the  rigid  body  motions  as  forced, 
sinusoidal,  damped  motions  - then  terms  appear  in  the 
differential  equations  which  can  produce  additional 
instabilities,  such  as  in  Mathieu's  equation. 

4.  In  the  case  of  general  motion  of  a flexible  body,  the 
combination  of  energy  methods  and  quasi  - coordinates 
can  produce  a practical  set  of  equations  that  govern  the 
aerodynamics,  flight  mechanics  and  structures  problems 
of  flexible  aircraft.  They  allow  the  determination  of  the 
effects  of  structural  flexibility  on  aircraft  performance, 
stability  and  air  loads  and  the  effects  of  the  "rigid  body" 
motions  on  aeroelastic  control  - effectiveness, 
divergence  and  flutter. 

5.  For  the  special  case  of  the  "flat  airplane"  with  perfect 
mass  distribution  and  perfect  modes,  a simple 
preliminary  estimate  of  the  effects  of  rigid  body 
motions  on  flexible  motions  would  be  to  replace  all  of 
the  structural  frequencies  by: 

t of-^wf-(pZ+Q 2) 
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Figure  1.  The  Sprung  Pendulum 
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Rgure  2.  Regions  of  Stability/Instability  for  Periodic 
Solutions  of  Mathieu's  Equation 
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Figure  3.  Instability  Region  for  Frequency  Ratios  near  1 
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Figure  10.  The  Effects  of  Acceleration  on  Hypersonic 
Divergence  and  Flutter 


